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摘要 ， 考虑 到 传统 Prony 算法 对 电力 系统 低频 振荡 模式 中 含有 衰减 较 快 的 分 量 分 
析 非 常 困 难以 及 容易 丢失 衰减 较 快 分 量 的 问题 ， 提 出 一 种 改进 的 快速 衰减 信号 的 提取 
方法 一 一 基于 Givens 变换 的 Prony 分 段 分 析 法 。 该 方法 将 采样 时 间 划 分 为 两 个 时 段 ， 
首先 分 析 第 二 时 段 求 慢 速 衰减 分 量 ， 然 后 分 析 第 一 时 段 求 快速 衰减 分 量 。 最 后 通过 算 
例 验证 了 该 算法 可 以 准确 提取 低频 振荡 模式 中 的 快速 衰减 分 量 ， 同 时 也 使 得 慢 速 衰减 
分 量 更 加 准确 。 
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Abstract: Considering that the traditional Prony algorithm is very difficult to analyze 


the fast attenuate component of the power system low frequency oscillation mode and the 
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Givens transform. This method divides the data acquisition time into two periods, and the 
second period is analyzed firstly to calculate the slow attenuate component, then the first 
period is analyzed for calculating the fast attenuate component. Numerical example proves 
that this algorithm can accurately extract fast attenuate component of low frequency 
oscillation mode, and also makes the slow attenuate component more accurate. 
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1 引言 


随 着 我 国 特 高 压 输电 网 络 的 引入 以 及 跨 大 区 域 
联网 进程 的 加 快 ， 低 频 振荡 "成 为 影响 电网 安全 稳 
定 的 重要 问题 ， 持 续 的 振荡 可 能 会 毁坏 电力 设备 继 
而 诱发 相继 事故 的 发 生 ， 造 成 严重 的 后 果 。 

实测 数据 分 析 表 明 ， 辨 识 振荡 模式 方法 有 短 时 傅 
里 叶 变 换 (Short-Time Fourier Transform，STEFT) 握 、 
小 波 分 析 法 0 、Prony 算法 和 HHT 变换 法 5 等 。 
综合 芳 虑 运算 量 、 耗 时 及 适用 性 ， 选 择 Prony 算法 ， 
该 算法 可 以 直接 从 时 域 中 通过 拟 合 求 得 信号 的 模 态 
信息 ， 不 需要 通过 频率 响应 求 取 ， 计 算 量 大 大 减少 ， 
适用 性 较 强 。 

为 了 更 有 效 地 分 析 系 统 的 低频 振荡 ， 提 取 全 部 
的 振东 特征， 本文 将 传统 的 Prony 算法 进行 改进 ， 
将 基于 Givens 变换 的 Prony 算法 进行 分 段 分 析 。 
Givens 变换 处 理 时 不 用 形成 法 方程 ， 可 以 直接 由 特 
征 方程 求 出 特征 方程 系数 ， 避 免 了 原 算法 求解 法 方 
程 困 难 和 误差 较 大 的 问题 ， 其 次 ,分 段 后 能 够 非常 
准确 地 辨识 出 低频 振东 模式 中 的 快速 衰减 分 量 ， 并 
且慢 速 训 减 分 量 也 更 加 精确 。 

2 Prony 算法 计算 步骤 

Prony 算法 是 提取 平稳 振荡 模式 的 常用 算 

法 。 它 针对 等 间距 采样 点 ， 假 设 模型 是 由 一 系列 具 


有 任意 振幅 、 相 位 、 频 率 和 衰减 因子 的 4 个 指数 函 
数 的 线性 组 合 ， 其 离散 时 间 的 函数 形式 为 


Xx(1) = When cos (2nft+0,) (1) 


式 中 ，4, 为 振幅 ， 0 为 初 相 (rad) ; fi 为 振荡 频率 
(Hz) ; Qi 为 衰减 因子 (阻尼 ) ; gq 为 信号 阶 数 。 

设 式 (1) 中 有 9 个 衰减 正 足 分 量 ， 然 后 将 余弦 
式 用 欧 拉 公式 表示 ， 即 


lr | 
cos (2nft+0,)= |e | (2) 


令 p=2g， 则 离散 时 间 的 函数 形式 可 表示 为 
Xm) -dbz n=0,1,…,N-1 (3) 


Prony 法 的 关键 是 要 认识 到 式 (1) 是 一 个 常 系 
数 线性 差分 方程 的 齐 次 解 。 即 有 


用 


ft(n)=- Sat(n-h) (4) 


i=] 


ChinaXiv 合 作 期 刊 


2018 年 12 月 


通过 一 系列 的 数学 变换 ， 可 得 到 差分 方程 式 为 


Rm) = -Yat(n-h) n=p,…,N-l (5) 
为 建立 Prony 算法 ， 定 义 实际 测量 值 x(n) 和 估 
计 值 #(n) 的 误差 为 e(n)， 即 


XxX(n)=X(n)+e(n) n=0,…,N-l1 (6) 


将 式 (5) 代入 式 (6)， 可 得 
Rm) = tln -Ah)+eln) 
-Saxn-b+ Yae(n-h) (7) 


由 式 (7) 可 以 看 出 ， 此 模型 是 一 个 由 误差 e(n) 
所 激励 的 自 回归 请 动 平 均 (Auto Regressive Moving 
Average) 模型 ARMA(p,p)， 它 具有 相同 的 自 回 
归 (Auto Regressive，AR) 和 滑动 平均 (Moving 
Average，MA) 参数 (a,) 。 为 使 模拟 信号 更 加 逼近 
于 真实 信号 ， 可 采用 平方 误差 最 小 原则 ， 即 关于 a， 


-au 4 的 最 小 二 乘 佑 计 使 得 了 |e(m)[ 值 最 小 ， 由 


此 将 得 到 一 组 有 一 定 求解 难度 的 非 线性 方程 。 
改变 误差 估计 的 目标 函数 可 以 避免 求解 非 线 性 


方程 ， 令 e(n)= > oueO- 昌 ， 则 式 (7) 可 变 为 
x(m) = ax(n-h)+eln) n=p,*…,N-1 (8) 


es N-l 
所 以 ， 变 史 |eoo 最 小 即 六 ecoj 最 小 ， 这 就 
n=p n=p 


是 扩展 的 Prony 算法 。 扩 展 的 Prony 算法 即 求解 下 
列 和 矩阵 方程 


xp) xp-D) x0) IT fs 
x(p+l) x(p) x) lol eptD 
x(N 1) x(N 2) 1 | 2 \e(N-1) 
(9) 
为 使 目标 函数 J(w = >|s(o)| 为 最 小 值 ， 令 
7(0) -0 ， 则 有 
90, 


$$ra-pro-nl-o i=1,…,p (10) 


n=p 
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对 应 的 最 小 误差 能 量 为 


p 
5,= Yo 


yx org (11) 


mp 
N-l 
定义 r(i,k)= Dx(n-k)x'(n-), i,j=0,…,p。 得 
n=p 
到 Prony 算法 的 法 方程 形式 为 


r(0,0) rz(0,1) :1 r(0,p) 1 2 


r(1,0) 7r(LD bpy 0 | i 


r(p,0) 7r(p,l) … r(p,p) N09,) \0 


求解 式 (12)， 可 得 参数 a1,…,a, 和 最 小 误差 能 
量 8, 的 估计 值 。 然 后 将 a1…,a, 代入 特征 方程 ， 进 
一 步 求 解 ， 得 到 特征 根 z,，i= 1,2,…,p。 则 式 (3) 
可 简化 为 关于 5b; 的 线性 方程 ， 用 和 矩阵 形式 表示 为 


Bb=X (13) 
其 中 
1 1 1 
| (14) 
el a 
Z| 22 2 
b=[b,…,b,] (15) 
则 
六 =[GaG] OB" K (17) 


式 中 ， 上 标 H 代表 共 轿 转 置 。 在 求解 过 程 中 ， 可 通 
过 下 述 关系 式 来 减少 计算 量 ， 提 高 计算 效率 ， 即 


Pu Pi en Bb, 
GG-=- 名 名 be (18) 
pb, Pb, pw 
其 中 
1 (19) 
(z12))-1 


式 (19) 中 ,“*” 代 表 共 力 。 
上 述 过 程 中 ， 由 于 事先 不 知道 低频 振荡 中 各 频 
率 分 量 的 个 数 ， 即 式 (9) 中 的 p 不 确定 ， 同 时 为 了 
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提高 算法 的 精度 ， 可 以 把 模型 的 阶 数 取 得 大 些 〈 即 
取 P。> p)， 再 通过 相应 方法 计算 出 方程 组 系数 矩阵 
的 有 效 秩 p。 

利用 式 (20) 计算 频率 f、 豪 减 因 子 a;/、 相 位 
9;、 振 幅 4,， 即 


f; =arctan [Im(z,)/ Re(z,)]/ (27At) 
o =In|z|/At 

0 =arctan [Im(b, )/ Re(2b, )] 

4 


(20) 


3 Prony 算法 的 改进 


3.1 传统 计算 方法 

在 利用 Prony 算法 对 电力 系统 低频 振荡 信号 进 
行 求解 的 过 程 中 ， 对 于 特征 方程 系数 的 求解 是 整个 
求解 过 程 的 核心 。 传 统 的 计算 方法 是 利用 Prony 
算法 的 法 方程 形式 进行 求解 ， 得 到 特征 方程 系数 a 
的 值 。 

将 式 (9) 简写 为 


Xa=d (21) 


式 中 , 怀 为 由 采样 数据 形成 的 系数 矩阵 ，a 为 特征 
方程 系数 ;4 为 信号 误差 。 
式 (21) 的 法 方程 形式 为 
X'Xa=X'd (22) 


利用 高 斯 消去 法 对 式 (22) 进行 求解 ， 求 得 特 
征 方程 系数 a。 

上 述 即 为 基于 法 方程 “1 形式 的 Prony 算法 的 
传统 计算 方法 。 
3.2 分 段 分 析 法 

现 选用 Givens 方法 对 特征 方程 系数 进行 求解 。 
利用 Givens 变换 王 对 方程 进行 求解 ， 得 到 特征 
方程 系数 时 ， 不 用 形成 法 方程 ， 可 以 直接 对 式 (21) 
进行 求解 。Givens 变换 方法 在 处 理 不 含 快速 衰减 分 
量 的 信号 时 ， 能 够 非常 准确 地 求 出 信号 的 参数 值 ， 
但 在 处 理 含 快速 衰减 分 量 的 信号 ， 尤 其 是 含有 两 个 
以 上 快速 衰减 分 量 的 信号 时 ，Givens 方法 及 其 他 传 
统 方法 不 一 定 能 够 取得 较 好 的 计算 结果 ， 有 时 甚至 
会 丢失 衰减 较 快 的 分 量 。 

为 解决 上 述 问题 ， 本 文 在 Givens 变换 方法 的 基 
础 上 ， 提 出 了 一 种 改进 方法 一 一 分 段 分 析 法 。 其 具 
体 实施 步骤 如 下 ; 
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(1) 读 入 采样 数据 。 设 采样 间隔 为 At， 采 样 时 
长 为 7。 为 了 能 够 精确 提取 快速 衰减 分 量 ， 要 求 采 
样 频率 在 0.01s 附近 取 值 ， 则 此 采样 数据 的 采样 间 
隔 为 At， 时 长 为 7， 采样 值 个 数 为 T/AL。 

(2) 时 段 划分 。 将 总 的 时 间 长 度 分 为 两 个 时 段 ， 
即 第 一 时 段 时 长 为 癌 ， 剩 余 时 段 划 为 第 二 时 段 ， 时 
长 为 7-t1。 第 一 时 段 和 第 二 时 段 的 分 割 点 为 A。 

(3) 第 二 时 段 Prony 法 低频 振荡 分 析 。 此 时 有 段 
分 析 慢 速 衰减 分 量 ， 采 样 间 隔 Ap 为 ,， 从 t, 时 刻 开 
始 每 间隔 从 采样 数据 中 取 一 个 采样 值 ， 即 每 隔 ty/ 
At 个 采样 数据 取 一 个 分 析 。 具 体 分 析 过 程 为 @@ 从 
t 时刻 开 始 每 间隔 从 采样 数据 中 取 一 个 采样 值 ， 
直到 了 结束 ，@ 对 经 过 消 品 处 理 后 的 采样 数据 进 
行 滤波 ， 滤 掉 信 号 中 的 高 频 成 分 ，@ 首 先 选取 阶 数 
q。，9。 大 于 有 效 秩 p， 然 后 根据 滤波 后 的 采样 数据 
形成 扩展 阶 矩 阵 及 .， 采 用 奇异 值 分 解法 对 扩展 阶 矩 
阵 R。 进行 计算 ， 求 出 扩展 阶 和 矩阵 的 所 有 奇异 值 ， 进 
而 确定 扩展 阶 和 矩阵 R. 的 有 效 秩 p; @ 重 新 选取 阶 数 
n，n 稍 大 于 有 效 秩 p， 小 于 q。.; @ 用 Givens 变换 求 
式 (21) 的 系数 a，@ 将 第 加 步 求 出 的 系数 代入 特 
征 方程 ， 求 出 特征 根 z,，@ 将 特征 根 代 入 式 (20)， 
求 出 频率 和 衰减 因子 ，@ 利 用 式 (17) ~ 式 (19) 
求 出 参数 b,，@ 将 b, 代 入 式 (20)， 求 出 幅 值 和 初 
相 ，@@ 最 后 将 四 个 参数 一 一 频率 、 训 减 因 子 、 幅 值 
和 初 相 代入 式 (1) 得 到 慢 速 衰减 分 量 信号 的 表达 式 
X(t) 

(4) 第 一 时 段 采样 数据 处 理 。 第 一 时 段 中 有 快 
速 衰减 分 量 和 慢 速 衰减 分 量 。 设 第 一 时 段 采样 间隔 
AN 为 4， 从 0s 开 始 进 行 采样 ， 直 到 + 时刻 为 止 。 
将 这 些 时 刻 点 代入 第 〈3) 步 已 经 求 出 的 慢 速 衰减 分 
量 表达 式 x(7?) 中 ， 得 到 这 一 时 间 段 每 个 时 刻 点 对 应 
的 慢 速 衰减 数据 ， 然 后 用 第 一 时 间 段 的 采样 数据 减 
去 刚才 求 出 的 慢 速 衰减 数据 ， 得 到 快速 衰减 数据 。 

(5) 第 一 时 段 Prony 法 低频 振荡 分 析 。 同 理 
采用 与 第 (3) 步 相 同 的 方法 分 析 已 经 得 到 的 快速 衰 
减 数据 ， 就 可 以 求 出 快速 衰减 分 量 的 四 个 模 态 信息 。 

下 图 为 本 文采 用 的 Prony 法 分 段 分 析 法 低频 振 
荡 模 式 分析 的 流程 图 。 

采用 分 段 分 析 法 后 ， 使 得 Givens 方法 在 对 含 快 
速 衰减 分 量 低频 信和 号 的 求解 效果 得 到 了 很 大 的 改善 。 
4 


算 例 分 析 


为 验证 本 文采 用 改进 的 Prony 法 分 段 分 析 法 能 
够 准确 地 提取 含 快速 衰减 分 量 的 低频 振荡 信号 ， 本 
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学 
读 入 采样 数据 
(经 过 消 噪 处 理 后 的 信号 ) 


Y 


滤波 ( 滤 掉 高 频 信 号 ) 


Y 


设 采样 时 间 间 隔 为 At， 采 样 时 长 为 7 


将 总 时 段 划 分 为 两 个 时 段 : 
第 一 时 段 为 含 快速 衰减 分 量 时 段 ， 第 二 时 段 为 慢 速 衰减 分 量 时 段 


第 二 时 段 Prony 法 低频 振 范 分析， 得 到 慢 速 衰减 分 量 信号 的 


达 式 x(?) 


Eee 

第 一 时 段 采 样 数据 处 理 ， 设 第 一 时 段 采 样 间隔 An 为 &， 从 0s 开 始 进行 

采样 至 4 时 刻 结束 。 将 这 些 时 刻 点 代入 z(D)， 得 到 该 时 间 段 每 个 时 刻 

点 对 应 的 慢 速 衰减 数据 ， 然 后 用 该 段 采样 数据 减 去 已 经 求 出 的 慢 速 
衰减 数据 ， 得 到 快速 衰减 数据 


第 一 时 段 采 用 Prony 法 低频 振荡 分 析 : 同 理 ， 采 用 与 第 二 时 段 相同 的 
分 析 方 法 便 可 得 到 信号 的 四 个 模 态 信息 


始 


Prony 法 分 段 分 析 法 低频 振荡 模式 分 析 流 程 图 
Fig. Flow chart of low frequency oscillation mode analysis 


by Prony method 


市 首先 将 本 改进 算法 和 基于 法 方程 形式 的 Prony 法 
分 段 分 析 算 法 进行 分 析 ， 其 次 将 本 改进 算法 与 基于 
Givens 变换 不 分 段 时 的 情况 进行 算 例 分 析 对 比 ， 最 
后 验证 了 本 改进 算法 具有 很 好 的 提取 效果 ， 能 够 十 
分 准确 地 提取 快速 衰减 分 量 。 

构造 输入 信号 为 


X(t)=10e "2 cos (2xx0.423t+110°)+ 
2e ”cos (2n x0.42t +13°)+ 
10e™™ cos (2x x0.2473t+20°)+ 
le ”cos (2nx1.0349t+60°)+ 
5.8e 8 cos (2a x 2.41 +10°)+ 
6.3e 5 cos (2n x 0.91+60°) (23) 


信号 模型 阶 数 n 统一 取 为 15。 上 述 信 号 包含 6 
个 正弦 分 量 ， 等 效 为 12 个 复 指数 分 量 。 对 上 述 信号 
每 0.01s 计算 一 个 值 ， 共 计算 1 000 个 值 ， 用 来 模拟 
采样 值 ， 则 此 采样 数据 的 采样 间隔 为 0.01s， 时 长 为 
10s， 采 样 值 个 数 为 1 000。 各 信和 号 分 量 参数 见 表 1。 
4.1 改进 算法 和 基于 法 方程 形式 的 Prony 分 段 分 析 

这 里 取信 号 分 段 点 为 1s。 第 一 时 段 分 析 快 速 
衰减 分 量 ， 采 样 时 间 间 隔 取 0.01s， 第 二 时 段 分 析 
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表 1 算 例 信号 的 参数 


Tab.1 The parameters of the example signal 


分 量 序号 ”频率 /Hz 豪 减 因 子 a 幅 值 4 初 相 9/(") 
1 0.4230 -0.0027 10.000 0 110.000 0 
2 0.4200 =0.265 2 2.000 0 13.000 0 
3 0.247 3 =0.031'1 10.0000 20.000 0 
4 1.034 9 -0.293 6 1.000 0 60.000 0 
Ss 2.400 0 -45.878 8 5.800 0 10.000 0 
0 0.900 0 =33. 则 36 6.300 0 60.000 0 


慢 速 衰减 分 量 ， 采 样 时 间 间 隔 取 0.1s。 表 2 为 基于 
Givens 变换 和 法 方程 形式 Prony 法 分 段 分 析 法 计算 
结果 比较 。 

表 2 两 种 算法 的 计算 结果 比较 
Tab.2 Comparison of calculation results of two kinds of 


algorithms 


分 量 序号 方法 频率 /Hz 衰减 因子 w 幅 值 4 初 相 9/(") 


Givens 变换 0.4230 -0.0027 10.0000 110.0000 


法 方程 形式 ”0.423 0 -0.0027 9.9999 109.9993 
Givens 变换 0.4200 -0.2652 2.0000 13.0000 
法 方程 形式 ”0.4200 =0.265 2 1.9998 13.0009 
Givens 变换 。 0.2473 -0.0311 10.0000 20.0000 
法 方程 形式 ”0.247 3 -0.0311 10.0000 20.0000 
Givens 变换 “1.0349 -0.293 6 1.0000 60.0000 
法 方程 形式 1.0349 -0.293 6 1.0000 60.0000 
Givens 变 换 2.4000 -45.8787 5.8001 9.999 1 
法 方程 形式 ”2.3659 -45.9230 6.2141 12.011:8 
Givens 变 换 0.9000 -55.1158 6.3001 60.0028 
° 法 方程 形式 ”1.5285 -56.5545 2.4217 60.1935 


由 表 2 可 见 ， 对 于 分 量 1 ~ 分 量 4 即 慢 速 衰减 分 
量 ， 基 于 Givens 变换 的 分 段 分 析 和 基于 法 方程 形式 
的 分 段 分 析 都 能 准确 地 提取 模 态 信息 ， 但 对 于 快速 衰 
减 分 量 ， 后 者 辨识 出 的 分 量 6 的 频率 为 1.528 5 Hz， 
而 理论 值 为 0.900 0Hz， 幅 值 为 2.421 7， 理 论 值 为 
6.300 1， 已 出 现 严 重 偏差 。 基 于 Givens 变换 的 分 段 
分 析 法 辨识 出 的 快速 衰减 分 量 精确 度 很 高 ， 与 理论 
值 误 差 很 小 。 

4.2 Givens 变换 方法 的 分 段 和 不 分 段 分 析 

信号 分 段 点 选 为 1s。 表 3 为 基于 Givens 变换 的 
分 段 和 不 分 段 分 析 的 计算 结果 比较 。 

由 表 3 可 见 ， 基 于 Givens 变换 的 分 段 法 能 十 分 
准确 地 提取 各 振荡 特征 的 低频 信号 的 模 态 信息 ， 不 
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表 3 Givens 变换 方法 的 分 段 和 不 分 段 计 算 结 果 比 较 
Tab.3 Comparison of the results of the segmentation and the 


whole calculation of the Givens transform method 


分 量 序号 ”方法 ”频率 fHz 衰减 因子 w 幅 值 4 和 初 相 0/(”) 


分 段 0.4230 -0.0027 10.0000 110.0000 
不 分 段 ”0.423 0 -0.0027 10.0000 110.0000 
分 段 0.4200 -0.2652 2.0000 13.0000 
? 不 分 段 0.4200 -0.2652 2.0000 ”13.0000 
分 段 ”0.2473 -0.0311 10.0000 20.0000 
: 不 分 段 ”0.2473 -0.0311 10.0000 20.0000 
分 段 1.0349 -0.2936 1.0000 60.0000 
“ 不 分 段 ” 1.0349 -0.2936 1.0000 60.0000 
分 段 2.4000 -45.8787 5.8001 9.999 1 
不 分 段 2.4266 -46.6631 7.6938 4.600 2 
6 分 段 0.9000 -55.1158 6.3001 60.002 8 


不 分 段 。 一 - - - 


分 段 时 ， 对 于 快速 衰减 分 量 6 却 无 法 辨识 。 该 算 例 
说 明 ， 本 文采 用 的 分 段 法 效果 较 好 ， 能 十 分 准确 地 
辨识 出 信号 中 的 快速 衰减 分 量 。 


5 结论 


(1) 在 振荡 模式 不 含 快速 衰减 分 量 情况 下 ， 传 
统 Prony 算法 能 够 较 准 确 地 辨识 振东 模 态 信息 。 当 
振 苞 模式 含有 两 个 及 以 上 快速 衰减 分 量 时 ， 传 统 
Prony 算法 由 于 采用 固定 的 采样 频率 ， 无 法 捕捉 豪 
减 较 快 的 分 量 ， 继 而 影响 慢 速 衰减 分 量 的 辨识 精度 。 

(2) 本 文 的 改进 算法 综合 了 Givens 变换 和 分 
段 Prony 法 的 优点 ， 对 处 理 消 噪 后 的 低频 振荡 信和 号 
有 独特 优势 。 采 用 Givens 变换 避免 了 传统 算法 求解 
法 方程 困难 的 问题 ， 提 高 了 数值 稳定 性 ， 减 少 了 计 
算 量 ， 利 用 分 段 分 析 设 置 不 同 的 分 辩 率 (采样 频率 ) 
避免 了 采用 同一 采样 频率 而 丢失 快速 衰减 分 量 。 该 
算法 占用 内 存 小 ， 能 非常 准确 地 提取 快速 衰减 分 量 ， 
适用 于 辨识 电力 系统 低频 振 范 模式 。 
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